This project seeks to estimate sport fish harvest and releases of rockfish in Alaska waters by improving on the Howard et al. (2020) methods and expand the time series back to 1977 when the statewide harvest survey (SWHS) was first implemented. This is essentially a Bayesian version of the Howard methods that allows for more appropriate and defensible sharing of information between areas, handles missing data in a more appropriate manor, accurately propagates uncertainty throughout the estimation procedure and replaces the Howard decision tree approach to low sample sizes with a hierarchical model. The methods and results for generating harvest estimates are generally consistent between the Bayesian model and the Howard methods. Harvest estimates are consistent with Howard estimates during contemporary times, but may differ based on more appropriate weighting of SWHS and logbook data, including estimating and correcting bias in the SWHS data.

Table 1. Summary of key improvements in reconstructiing sport fish removals of rockfish using the Bayesian model as compared to the Howard et al. (2020) methods.
Issue Howard Bayes
Time series 1999 - present 1977 - present
Bias in SWHS Not explicitly dealt with. Relies on logbook data and ratios of guided/unguided from SWHS data to estimate unguided releases and harvests. Explicitly estimates bias in SWHS harvest and release estimates based on logbook data.
Species composition of releases Assumes that species composition of releases is equal to that of the harvest, which is not evident in the logbook data. Recognizes different release probabilities by species / species assemblage and estimates it from logbook data and bias corrected SWHS data
Sample size limitations Uses sample size threshholds such that when areas fall below those threshholds values are borrowed from nearby areas. Uses a hierarchichacal modelling approach that shares information between areas in the same region. Thus all data is used, even with small sample sizes. This is a more sound method that avoids assumptions and uses all of the data.
Error propogation Error is propogated when variance estimates are available, but there is uncertainty associated with borrowing values from nearby areas, or the assumption of species compositions being identical in harvest and releases, are not dealt with. By breaking the assumption that species composition is equal between harvests and releases, uncertainty in the release estimates is more reflective of the fishery. Furthermore, the hyerarchichal approach more accurately captures uncertainy within and between areas within a region.
**Figure 1.**- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.

Figure 1.- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.


**Figure 2.**- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).

Figure 2.- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).


Results

**Figure X.**- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Figure X.- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Estimate comparison

Since previous estimates of rockfish harvest have been produced these first 3 graphs will be used to show how the modeled estimates compare to the estimates produced earlier. For total rockfish the estimates are in general agreement although differences are noted. These estimates should be more reliable because they include both SWHS and guide logbook data, handle variance more appropriately, use hierarchical distributions when data is missing, directly consider observation error and are produced using reproducible research.

**Figure 3.**- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 3.**- Total rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


Notes from Adam: When looking at only black rockfish the most significant differences are for the Prince William Sound Inside area. I did not spend a great deal of time tracking this down although it looks like the previous version used bad values for \(P_{(black)ayu}\) for at least unguided anglers. For the moment I would ignore the results for BSIA and SOKO2SAP. I think it is possible to give approximate values for these areas but it will require a little more coding which I have yet to do.

**Figure 4.**- Black rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 4.- Black rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


And black rockfish releases…

**Figure 5.**- Black rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 5.- Black rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 6.**- Yellow rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 6.- Yellow rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 7.**- Yellow rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 7.- Yellow rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 8.**- DSR rockfish (excluding yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 8.- DSR rockfish (excluding yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 9.**- DSR rockfish releases (including yelloweye) 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 9.- DSR rockfish releases (including yelloweye) 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 11.**- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 11.- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 12.**- Slope rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 12.- Slope rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

**Figure 13.**- *Other* rockfish (non-pelagic and non-yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 13.- Other rockfish (non-pelagic and non-yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 14.**- *Other* rockfish (non-pelagic and non-yelloweye) releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 14.- Other rockfish (non-pelagic and non-yelloweye) releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



Total Biomass Removal Estimates

**Figure 15.**- Black rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 15.- Black rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.



**Figure 16.**- Yellow rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 16.- Yellow rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

**Figure 17.**- Pelagic rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 17.- Pelagic rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 18.**- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 18.- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 19.**- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 19.- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


Model fit

Logbook residuals

**Figure 20.**- Residuals from logbook harvests.

Figure 20.- Residuals from logbook harvests.


SWHS residuals

**Figure 21.**- Residuals from SWHS harvests.

Figure 21.- Residuals from SWHS harvests.



**Figure 22.**- Residual of SWHS releases.

Figure 22.- Residual of SWHS releases.

Parameter estimates

P(Charter)

These histograms show the posterior distribution of the mean percent of rockfish harvested by the charter fleet.

**Figure 23.**- Mean percent of harvest by charter anglers.

Figure 23.- Mean percent of harvest by charter anglers.


When considered annually we see the percent of rockfish harvested by the charter fleet follows our data fairly well although the model smooths out the changes and we just do not have much information about this ratio. Prior to 2011 the percent charter is confounded with SWHS bias and should be mostly discounted.

**Figure 24.**- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

Figure 24.- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

P(Harvest)

These plots show the fitted logistic line to the proportion of caught rockfish that are harvested. These estimates are used for hindcasting catch estimates based on the harvest data in early years when catch estimates are unavailable.


**Figure 25.**- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 25.- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 26.**- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 26.- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 27.**- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.

Figure 27.- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.


**Figure 28.**- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 28.- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.


**Figure 29.**- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 29.- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

SWHS bias

Figure 23 shows the mean estimate for SWHS bias in harvests and releases. Cook Inlet, North Gulf Coast and North Southeast Inside all look pretty good while most other areas have substantial bias. Prince William Sound Inside has the largest bias. Bias in release estimates is substantial and whereas the SWHS appears to underestimate harvests, it appears to greatly overestimates releases by a factor of 2 or more in most areas as derived from logbook reported releases.

**Figure 30.**- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 30.- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.


Our estimates of SWHS harvest bias track observations fairly well when he have guided harvest estimates. The estimates of release bias in the SWHS data track observed patterns to an extent, but appear to smooth these more volatile disagreements with the logbook data. Adam postulated in his initial start on this that some of this could be the result of the estimates of the proportion guided. This value was not modelled with a trend and thus applies a constant estimate when hindcasting. Data on these relationships could greatly improve this model.

**Figure 31.**- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 31.- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

P(pelagic)

We model the percentage of pelagic rockfish in the harvest because we have the information for charter anglers (via logbooks) starting in 1998. Other than looking at the model estimates you can use Figure 25 to compare the two data streams for pelagic rockfish harvest. In general they are in agreement with major exceptions in Price William Sound inside, Prince William Sound outside (early in the time series) and South Southeast inside.

**Figure 32.**- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

Figure 32.- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

P(black|pelagic)

Note that in Southeast Alaska we only have composition data starting in 2006. Tania dug up old SE data, but it did not provide any useful data for species apportionment. For the most part, P(black|pelagic) is relatively constant across areas, with the exception of Cook Inlet and NSEI in Southeast AK. It may be worth discussing whether the shifts in those areas is a result of improved or changing species identification rather than actual shift in the species composition of the catch.

**Figure 33.**- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

Figure 33.- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

P(yelloweye|non-pelagic / yelloweye|DSR)

**Figure 34.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 34.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(DSR|non-pelagic)

**Figure 35.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

Figure 35.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

P(slope|non-pelagic)

**Figure 36.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 36.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.



P(slope|non-pelagic & non-yellowye) For release estimates

## NULL



Weight Fits

**Figure 37.**- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 37.- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 38.**- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 38.- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 39.**- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 39.- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 40.**- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 40.- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 41.**- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 41.- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


Summary of unconverged parameters:

The proportion of unconverged species composition proportional estimates that are near the bounds of 0 or 1 is 0.4612245

The proportion of unconverged true parameters (defined as Rhat > 1.01) that are also in excess of an Rhat of 1.1 is 0.08 while 0.92 are less than 1.1.

Table Xa. Species proportion derived quantities that are not with 0.05 of the 0 and 1 bounds and are not converged relative to an Rhat of 1.01.

variable

area_n

year

user

mean

median

lower

upper

Rhat

sd

skew

area

p_slope

11

22

2

0.098

0.074

0.013

0.326

1.013

0.085

0.283

CSEO

p_black

11

23

2

0.868

0.904

0.572

0.985

1.012

0.112

-0.322

CSEO

p_pelagic

13

1

1

0.269

0.249

0.064

0.570

1.040

0.134

0.149

NSEI

p_pelagic

13

1

2

0.110

0.093

0.020

0.298

1.036

0.074

0.230

NSEI

p_yellow

13

1

1

0.394

0.356

0.128

0.926

1.048

0.191

0.199

NSEI

p_yellow

13

1

2

0.397

0.359

0.127

0.929

1.053

0.193

0.197

NSEI

p_pelagic

13

10

1

0.294

0.278

0.081

0.591

1.021

0.134

0.119

NSEI

p_pelagic

13

10

2

0.120

0.107

0.028

0.293

1.030

0.069

0.188

NSEI

p_yellow

13

10

1

0.396

0.356

0.126

0.923

1.046

0.194

0.206

NSEI

p_yellow

13

10

2

0.398

0.355

0.126

0.929

1.046

0.195

0.221

NSEI

p_pelagic

13

11

1

0.294

0.278

0.090

0.570

1.023

0.129

0.124

NSEI

p_pelagic

13

11

2

0.122

0.107

0.029

0.296

1.021

0.070

0.216

NSEI

p_yellow

13

11

1

0.396

0.357

0.133

0.926

1.041

0.193

0.202

NSEI

p_yellow

13

11

2

0.394

0.356

0.132

0.926

1.039

0.195

0.195

NSEI

p_pelagic

13

12

1

0.305

0.289

0.097

0.591

1.012

0.130

0.123

NSEI

p_pelagic

13

12

2

0.125

0.111

0.031

0.298

1.024

0.071

0.197

NSEI

p_yellow

13

12

1

0.395

0.353

0.126

0.934

1.044

0.195

0.216

NSEI

p_yellow

13

12

2

0.394

0.356

0.127

0.930

1.040

0.193

0.196

NSEI

p_pelagic

13

13

1

0.309

0.296

0.099

0.594

1.021

0.129

0.101

NSEI

p_pelagic

13

13

2

0.129

0.113

0.031

0.312

1.020

0.073

0.220

NSEI

p_yellow

13

13

1

0.393

0.352

0.127

0.932

1.036

0.194

0.211

NSEI

p_yellow

13

13

2

0.395

0.357

0.123

0.932

1.040

0.194

0.196

NSEI

p_pelagic

13

14

1

0.301

0.290

0.107

0.561

1.013

0.120

0.092

NSEI

p_pelagic

13

14

2

0.129

0.115

0.035

0.302

1.022

0.070

0.201

NSEI

p_yellow

13

14

1

0.395

0.361

0.128

0.920

1.040

0.193

0.176

NSEI

p_yellow

13

14

2

0.397

0.356

0.133

0.925

1.036

0.194

0.211

NSEI

p_pelagic

13

15

1

0.274

0.258

0.092

0.522

1.018

0.114

0.141

NSEI

p_pelagic

13

15

2

0.126

0.112

0.035

0.297

1.024

0.068

0.207

NSEI

p_yellow

13

15

1

0.390

0.351

0.125

0.931

1.038

0.193

0.202

NSEI

p_yellow

13

15

2

0.392

0.350

0.131

0.919

1.045

0.194

0.216

NSEI

p_pelagic

13

16

2

0.137

0.123

0.039

0.310

1.015

0.071

0.198

NSEI

p_yellow

13

16

1

0.391

0.353

0.126

0.929

1.035

0.194

0.195

NSEI

p_yellow

13

16

2

0.398

0.359

0.127

0.925

1.038

0.193

0.202

NSEI

p_pelagic

13

17

2

0.142

0.127

0.041

0.319

1.014

0.072

0.207

NSEI

p_yellow

13

17

1

0.394

0.357

0.131

0.933

1.046

0.193

0.192

NSEI

p_yellow

13

17

2

0.393

0.353

0.130

0.927

1.043

0.192

0.208

NSEI

p_yellow

13

18

1

0.395

0.358

0.124

0.931

1.036

0.195

0.190

NSEI

p_yellow

13

18

2

0.397

0.359

0.131

0.929

1.043

0.193

0.197

NSEI

p_pelagic

13

19

1

0.338

0.330

0.140

0.591

1.010

0.117

0.069

NSEI

p_yellow

13

19

1

0.392

0.354

0.124

0.927

1.042

0.195

0.195

NSEI

p_yellow

13

19

2

0.396

0.356

0.130

0.923

1.044

0.194

0.206

NSEI

p_pelagic

13

2

1

0.269

0.250

0.065

0.573

1.033

0.134

0.141

NSEI

p_pelagic

13

2

2

0.109

0.093

0.022

0.288

1.037

0.071

0.226

NSEI

p_yellow

13

2

1

0.390

0.346

0.130

0.939

1.042

0.195

0.226

NSEI

p_yellow

13

2

2

0.398

0.358

0.127

0.926

1.038

0.194

0.206

NSEI

p_yellow

13

20

1

0.391

0.352

0.126

0.930

1.042

0.193

0.203

NSEI

p_yellow

13

20

2

0.395

0.358

0.130

0.924

1.047

0.191

0.193

NSEI

p_yellow

13

21

1

0.394

0.358

0.127

0.928

1.044

0.191

0.188

NSEI

p_yellow

13

21

2

0.396

0.353

0.122

0.934

1.042

0.195

0.221

NSEI

p_yellow

13

22

1

0.393

0.354

0.127

0.922

1.043

0.194

0.201

NSEI

p_yellow

13

22

2

0.397

0.358

0.133

0.927

1.043

0.193

0.202

NSEI

p_yellow

13

23

1

0.356

0.313

0.098

0.915

1.047

0.194

0.221

NSEI

p_yellow

13

23

2

0.446

0.417

0.150

0.939

1.038

0.192

0.151

NSEI

p_yellow

13

24

1

0.361

0.321

0.096

0.905

1.039

0.197

0.203

NSEI

p_yellow

13

24

2

0.447

0.418

0.155

0.934

1.032

0.192

0.151

NSEI

p_yellow

13

25

1

0.363

0.320

0.106

0.906

1.044

0.193

0.223

NSEI

p_yellow

13

25

2

0.441

0.413

0.155

0.926

1.027

0.192

0.146

NSEI

p_yellow

13

26

1

0.353

0.311

0.102

0.898

1.045

0.191

0.220

NSEI

p_yellow

13

26

2

0.439

0.408

0.155

0.933

1.030

0.191

0.162

NSEI

p_yellow

13

27

1

0.346

0.305

0.096

0.886

1.044

0.190

0.216

NSEI

p_yellow

13

27

2

0.436

0.404

0.151

0.924

1.028

0.191

0.168

NSEI

p_yellow

13

28

1

0.354

0.315

0.103

0.865

1.047

0.184

0.212

NSEI

p_yellow

13

28

2

0.438

0.413

0.150

0.897

1.028

0.187

0.133

NSEI

p_yellow

13

29

1

0.341

0.304

0.105

0.826

1.036

0.178

0.208

NSEI

p_yellow

13

29

2

0.434

0.411

0.153

0.875

1.027

0.179

0.128

NSEI

p_pelagic

13

3

1

0.272

0.254

0.067

0.568

1.040

0.135

0.134

NSEI

p_pelagic

13

3

2

0.110

0.091

0.022

0.299

1.031

0.072

0.265

NSEI

p_yellow

13

3

1

0.395

0.356

0.126

0.926

1.042

0.195

0.200

NSEI

p_yellow

13

3

2

0.393

0.352

0.129

0.930

1.053

0.195

0.210

NSEI

p_pelagic

13

4

1

0.273

0.254

0.066

0.571

1.040

0.133

0.143

NSEI

p_pelagic

13

4

2

0.110

0.094

0.023

0.287

1.030

0.071

0.226

NSEI

p_yellow

13

4

1

0.395

0.354

0.137

0.925

1.052

0.194

0.211

NSEI

p_yellow

13

4

2

0.400

0.365

0.131

0.922

1.045

0.192

0.182

NSEI

p_pelagic

13

5

1

0.278

0.261

0.074

0.569

1.027

0.133

0.128

NSEI

p_pelagic

13

5

2

0.112

0.095

0.023

0.292

1.044

0.070

0.243

NSEI

p_yellow

13

5

1

0.396

0.358

0.129

0.922

1.037

0.195

0.195

NSEI

p_yellow

13

5

2

0.396

0.362

0.128

0.927

1.042

0.192

0.177

NSEI

p_pelagic

13

6

1

0.279

0.261

0.075

0.594

1.036

0.134

0.134

NSEI

p_pelagic

13

6

2

0.112

0.097

0.024

0.291

1.045

0.070

0.214

NSEI

p_yellow

13

6

1

0.390

0.349

0.126

0.937

1.042

0.193

0.212

NSEI

p_yellow

13

6

2

0.401

0.363

0.126

0.934

1.045

0.196

0.194

NSEI

p_pelagic

13

7

1

0.281

0.262

0.071

0.577

1.032

0.132

0.143

NSEI

p_pelagic

13

7

2

0.114

0.099

0.026

0.290

1.030

0.069

0.218

NSEI

p_yellow

13

7

1

0.391

0.350

0.132

0.921

1.042

0.193

0.213

NSEI

p_yellow

13

7

2

0.397

0.361

0.125

0.928

1.041

0.193

0.187

NSEI

p_pelagic

13

8

1

0.287

0.271

0.080

0.586

1.027

0.134

0.119

NSEI

p_pelagic

13

8

2

0.115

0.099

0.027

0.294

1.029

0.070

0.228

NSEI

p_yellow

13

8

1

0.397

0.362

0.124

0.926

1.039

0.196

0.179

NSEI

p_yellow

13

8

2

0.399

0.362

0.129

0.922

1.039

0.195

0.190

NSEI

p_pelagic

13

9

1

0.290

0.275

0.084

0.567

1.032

0.130

0.115

NSEI

p_pelagic

13

9

2

0.118

0.103

0.027

0.298

1.034

0.071

0.213

NSEI

p_yellow

13

9

1

0.396

0.355

0.125

0.931

1.040

0.196

0.210

NSEI

p_yellow

13

9

2

0.397

0.360

0.131

0.926

1.041

0.194

0.191

NSEI

p_pelagic

16

1

1

0.370

0.359

0.133

0.652

1.031

0.138

0.080

SSEO

p_pelagic

16

1

2

0.347

0.337

0.128

0.628

1.023

0.132

0.076

SSEO

p_pelagic

16

10

1

0.374

0.365

0.141

0.657

1.027

0.134

0.067

SSEO

p_pelagic

16

10

2

0.352

0.337

0.127

0.628

1.020

0.133

0.113

SSEO

p_pelagic

16

11

1

0.372

0.366

0.146

0.638

1.024

0.132

0.045

SSEO

p_pelagic

16

11

2

0.354

0.344

0.139

0.639

1.020

0.132

0.076

SSEO

p_pelagic

16

12

1

0.377

0.369

0.146

0.662

1.022

0.135

0.059

SSEO

p_pelagic

16

12

2

0.354

0.344

0.133

0.628

1.022

0.132

0.076

SSEO

p_pelagic

16

13

1

0.377

0.368

0.155

0.652

1.026

0.130

0.069

SSEO

p_pelagic

16

13

2

0.356

0.350

0.134

0.625

1.018

0.130

0.046

SSEO

p_pelagic

16

14

1

0.366

0.356

0.145

0.642

1.017

0.131

0.076

SSEO

p_pelagic

16

14

2

0.356

0.348

0.141

0.624

1.022

0.128

0.063

SSEO

p_pelagic

16

15

1

0.369

0.360

0.152

0.638

1.017

0.125

0.072

SSEO

p_pelagic

16

15

2

0.360

0.350

0.150

0.630

1.016

0.125

0.080

SSEO

p_pelagic

16

16

1

0.363

0.356

0.152

0.627

1.018

0.126

0.056

SSEO

p_pelagic

16

16

2

0.362

0.351

0.146

0.624

1.014

0.127

0.087

SSEO

p_pelagic

16

17

1

0.424

0.418

0.192

0.677

1.013

0.127

0.047

SSEO

p_pelagic

16

17

2

0.380

0.369

0.156

0.650

1.014

0.130

0.085

SSEO

p_pelagic

16

18

1

0.403

0.397

0.179

0.662

1.016

0.125

0.048

SSEO

p_pelagic

16

18

2

0.379

0.369

0.155

0.648

1.013

0.130

0.077

SSEO

p_pelagic

16

19

1

0.394

0.385

0.174

0.655

1.015

0.125

0.072

SSEO

p_pelagic

16

2

1

0.368

0.357

0.139

0.655

1.022

0.138

0.080

SSEO

p_pelagic

16

2

2

0.349

0.336

0.125

0.630

1.017

0.133

0.098

SSEO

p_pelagic

16

21

2

0.380

0.375

0.170

0.627

1.012

0.121

0.041

SSEO

p_pelagic

16

3

1

0.371

0.364

0.138

0.647

1.028

0.136

0.051

SSEO

p_pelagic

16

3

2

0.348

0.340

0.122

0.626

1.026

0.133

0.060

SSEO

p_pelagic

16

4

1

0.367

0.353

0.133

0.658

1.020

0.138

0.102

SSEO

p_pelagic

16

4

2

0.352

0.342

0.129

0.631

1.017

0.132

0.076

SSEO

p_pelagic

16

5

1

0.372

0.363

0.136

0.659

1.018

0.135

0.067

SSEO

p_pelagic

16

5

2

0.351

0.343

0.130

0.626

1.024

0.131

0.061

SSEO

p_pelagic

16

6

1

0.372

0.365

0.133

0.655

1.022

0.136

0.051

SSEO

p_pelagic

16

6

2

0.348

0.336

0.131

0.628

1.021

0.132

0.091

SSEO

p_pelagic

16

7

1

0.370

0.360

0.139

0.645

1.023

0.135

0.074

SSEO

p_pelagic

16

7

2

0.351

0.340

0.129

0.628

1.022

0.133

0.083

SSEO

p_pelagic

16

8

1

0.371

0.362

0.145

0.652

1.019

0.135

0.067

SSEO

p_pelagic

16

8

2

0.353

0.340

0.136

0.631

1.025

0.132

0.098

SSEO

p_pelagic

16

9

1

0.372

0.362

0.143

0.662

1.021

0.135

0.074

SSEO

p_pelagic

16

9

2

0.353

0.344

0.133

0.628

1.019

0.130

0.069

SSEO

p_black

8

40

1

0.461

0.413

0.345

0.960

1.168

0.159

0.302

afognak

Table Xb. True parameters with an Rhat in excess of 1.1

variable

area_n

year/species

mean

median

lower

upper

Rhat

sd

skew

year

area

beta2_pH

10

3

0.302

0.225

0.099

0.931

1.202

0.352

0.219

northeast

beta0_pH

11

1

0.136

0.140

-0.409

0.676

1.101

0.280

-0.014

CSEO

beta2_pH

11

1

0.477

0.391

0.206

1.237

1.129

0.306

0.281

CSEO

beta0_pelagic

13

-1.246

-1.179

-2.637

-0.202

1.108

0.633

-0.106

NSEI

beta1_pelagic

13

3.399

3.406

1.738

4.884

1.108

0.802

-0.009

NSEI

beta1_yellow

13

0.917

0.717

0.023

3.865

1.142

0.892

0.224

NSEI

beta0_pH

15

4

0.708

0.813

-0.810

1.599

1.228

0.687

-0.153

SSEI

beta1_pH

15

4

1.063

0.977

0.115

2.645

1.203

0.754

0.114

SSEI

tau_beta2_pH

2

3

0.269

0.112

0.041

1.226

1.173

1.010

0.155

NG

tau_beta1_pH

3

4

4.207

0.441

0.011

34.738

1.124

10.280

0.366

PWSI

beta2_pH

8

1

0.573

0.278

0.104

3.388

1.231

1.067

0.276

afognak

beta0_H

8

6.501

6.277

-20.023

29.978

1.135

38.526

0.006

afognak

Table Xc. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.01.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.01

405.00

205.00

925.00

324.00

Prop. of total

0.66

0.34

0.74

0.26

Prop. median less than 1

0.29

0.22

0.34

0.42

Prop. median less than 10

0.38

0.54

0.54

0.72

Prop. median less than 50

0.52

0.79

0.67

0.78

Median of the median values

41.45

7.74

6.40

1.80

Mean of the median values

440.82

277.07

298.53

140.67

Mean skew of posteriors

0.25

0.29

0.27

0.24

Table Xd. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.1.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.1

1.00

5.00

48.00

36.00

Prop. of total

0.17

0.83

0.57

0.43

Prop. median less than 1

0.00

0.40

0.77

0.58

Prop. median less than 10

0.00

0.80

0.85

0.89

Prop. median less than 50

0.00

1.00

0.96

0.89

Median of the median values

1,866.71

2.41

0.01

0.25

Mean of the median values

1,866.71

9.87

13.21

13.64

Mean skew of posteriors

0.30

0.20

0.14

0.18

Table Xe. All parameters and derived quantities that are not converged relative to an Rhat of 1.1.

variable

Freq

beta0_H

1

beta0_pH

1

beta1_pH

1

beta1_yellow

1

beta2_pH

3

eps_bl

1

H_ay

1

H_ayu

1

Hb_ay

2

Hb_ayg

1

Hb_ayu

1

Ho_ayg

1

Hp_ay

1

Hp_ayu

1

Hs_ayu

3

Hy_ayg

2

p_black

1

p_slope

1

R_ay

1

R_ayg

2

R_ayu

4

Rb_ay

1

Rb_ayg

3

Rb_ayu

4

Rd_ay

4

Rd_ayg

6

Rd_ayu

3

re_black

1

re_pH

1

Ro_ay

5

Ro_ayg

7

Ro_ayu

26

Rp_ay

1

Rp_ayg

2

Rp_ayu

4

Rs_ayu

8

Ry_ay

12

Ry_ayg

4

Ry_ayu

23

tau_beta1_pH

1

tau_beta2_pH

1

Parameter estimates:

Summary Table of Parameter Estimates
Parameter mean sd Lower_CI Median Upper_CI
mu_bc_H[1] -0.016 0.070 -0.147 -0.020 0.125
mu_bc_H[2] 0.074 0.054 -0.025 0.073 0.183
mu_bc_H[3] -0.269 0.067 -0.394 -0.271 -0.134
mu_bc_H[4] -0.929 0.184 -1.305 -0.926 -0.578
mu_bc_H[5] 1.323 1.141 0.067 1.053 4.172
mu_bc_H[6] -2.401 0.317 -3.026 -2.406 -1.776
mu_bc_H[7] -0.452 0.099 -0.651 -0.449 -0.264
mu_bc_H[8] 0.258 0.366 -0.351 0.220 1.084
mu_bc_H[9] -0.248 0.123 -0.490 -0.250 -0.004
mu_bc_H[10] 0.059 0.065 -0.060 0.055 0.194
mu_bc_H[11] -0.125 0.037 -0.199 -0.125 -0.056
mu_bc_H[12] -0.251 0.107 -0.476 -0.246 -0.045
mu_bc_H[13] -0.152 0.072 -0.287 -0.153 -0.003
mu_bc_H[14] -0.299 0.090 -0.474 -0.296 -0.126
mu_bc_H[15] -0.326 0.046 -0.415 -0.326 -0.232
mu_bc_H[16] -0.167 0.374 -0.793 -0.198 0.625
mu_bc_R[1] 1.042 0.121 0.798 1.044 1.271
mu_bc_R[2] 1.011 0.078 0.847 1.012 1.160
mu_bc_R[3] 0.225 0.089 0.044 0.227 0.404
mu_bc_R[4] 0.655 0.199 0.245 0.663 1.028
mu_bc_R[5] 1.666 0.363 0.977 1.658 2.404
mu_bc_R[6] -1.436 0.338 -2.087 -1.448 -0.759
mu_bc_R[7] 0.517 0.167 0.180 0.517 0.840
mu_bc_R[8] 0.640 0.146 0.361 0.639 0.927
mu_bc_R[9] 0.378 0.168 0.034 0.383 0.682
mu_bc_R[10] 0.895 0.101 0.695 0.896 1.098
mu_bc_R[11] 1.110 0.075 0.963 1.111 1.260
mu_bc_R[12] 0.892 0.163 0.566 0.894 1.220
mu_bc_R[13] 1.060 0.094 0.871 1.061 1.242
mu_bc_R[14] 0.865 0.120 0.620 0.864 1.091
mu_bc_R[15] 0.732 0.094 0.543 0.733 0.914
mu_bc_R[16] 1.264 0.115 1.035 1.265 1.489
tau_pH[1] 2.454 0.189 2.111 2.446 2.849
tau_pH[2] 1.582 0.149 1.307 1.577 1.881
tau_pH[3] 2.268 0.300 1.720 2.256 2.897
tau_pH[4] 1.048 0.142 0.792 1.041 1.348
tau_pH[5] 1.015 0.153 0.745 1.002 1.351
beta0_pH[1,1] 0.567 0.179 0.206 0.571 0.895
beta0_pH[2,1] 1.426 0.163 1.080 1.434 1.725
beta0_pH[3,1] 1.401 0.222 0.880 1.427 1.763
beta0_pH[4,1] 1.560 0.225 1.087 1.567 1.984
beta0_pH[5,1] -0.834 0.380 -1.744 -0.797 -0.205
beta0_pH[6,1] -0.541 0.435 -1.669 -0.437 0.032
beta0_pH[7,1] 0.752 0.011 0.732 0.752 0.774
beta0_pH[8,1] -0.710 0.445 -1.920 -0.645 -0.004
beta0_pH[9,1] -0.344 0.305 -1.071 -0.311 0.155
beta0_pH[10,1] 0.264 0.187 -0.156 0.279 0.596
beta0_pH[11,1] 0.136 0.280 -0.409 0.140 0.676
beta0_pH[12,1] 0.646 0.195 0.240 0.651 1.004
beta0_pH[13,1] -0.655 0.149 -0.967 -0.647 -0.376
beta0_pH[14,1] 0.015 0.171 -0.333 0.018 0.343
beta0_pH[15,1] 0.048 0.198 -0.381 0.067 0.390
beta0_pH[16,1] -0.866 0.202 -1.294 -0.855 -0.501
beta0_pH[1,2] 2.567 0.082 2.410 2.566 2.729
beta0_pH[2,2] 3.079 0.022 3.036 3.079 3.122
beta0_pH[3,2] 3.626 0.041 3.548 3.627 3.704
beta0_pH[4,2] 2.936 0.019 2.899 2.936 2.974
beta0_pH[5,2] 4.555 1.509 2.774 4.267 7.977
beta0_pH[6,2] 3.324 0.158 3.028 3.320 3.644
beta0_pH[7,2] 2.023 0.049 1.931 2.022 2.123
beta0_pH[8,2] 2.824 0.066 2.698 2.823 2.955
beta0_pH[9,2] 3.961 0.208 3.569 3.948 4.387
beta0_pH[10,2] 4.002 0.129 3.753 4.002 4.254
beta0_pH[11,2] 1.888 0.009 1.871 1.888 1.906
beta0_pH[12,2] 1.803 0.068 1.676 1.802 1.940
beta0_pH[13,2] 2.498 0.026 2.446 2.498 2.549
beta0_pH[14,2] 1.805 0.025 1.757 1.804 1.854
beta0_pH[15,2] 2.550 0.017 2.518 2.550 2.584
beta0_pH[16,2] 2.776 0.025 2.728 2.776 2.826
beta0_pH[1,3] 0.686 0.614 -1.059 0.876 1.372
beta0_pH[2,3] 2.126 0.018 2.093 2.126 2.160
beta0_pH[3,3] 2.690 0.024 2.644 2.689 2.738
beta0_pH[4,3] 2.981 0.034 2.914 2.980 3.048
beta0_pH[5,3] 3.745 2.598 0.873 3.065 10.672
beta0_pH[6,3] 1.255 0.327 0.686 1.239 1.952
beta0_pH[7,3] 0.856 0.049 0.761 0.856 0.952
beta0_pH[8,3] 0.287 0.120 0.026 0.294 0.507
beta0_pH[9,3] -2.863 1.118 -5.347 -2.717 -0.977
beta0_pH[10,3] -0.913 0.906 -2.976 -0.783 0.568
beta0_pH[11,3] -0.194 11.722 -23.349 -0.226 23.099
beta0_pH[12,3] 0.025 11.859 -22.477 0.104 23.377
beta0_pH[13,3] -0.027 11.604 -22.397 0.112 23.104
beta0_pH[14,3] 0.001 11.692 -22.808 0.065 23.858
beta0_pH[15,3] -0.098 11.782 -22.650 0.031 23.085
beta0_pH[16,3] -0.023 11.621 -22.980 -0.100 23.677
beta0_pH[11,4] 1.594 0.636 0.236 1.626 2.816
beta0_pH[12,4] 1.505 1.077 -0.987 1.500 3.559
beta0_pH[13,4] 1.983 0.804 0.114 2.075 3.290
beta0_pH[14,4] 0.708 0.613 -0.697 0.759 1.737
beta0_pH[15,4] 0.708 0.687 -0.810 0.813 1.599
beta0_pH[16,4] 1.308 0.815 -0.540 1.384 2.732
beta0_pH[11,5] 1.287 0.105 1.104 1.280 1.511
beta0_pH[12,5] -2.302 0.428 -3.232 -2.260 -1.582
beta0_pH[13,5] 0.761 0.055 0.660 0.757 0.878
beta0_pH[14,5] 0.643 0.237 0.209 0.633 1.139
beta0_pH[15,5] 2.200 0.417 1.576 2.129 3.269
beta0_pH[16,5] -0.769 0.068 -0.895 -0.772 -0.632
beta1_pH[1,1] 2.543 0.493 1.745 2.476 3.709
beta1_pH[2,1] 1.562 0.304 0.999 1.558 2.194
beta1_pH[3,1] 1.800 0.442 1.103 1.735 2.944
beta1_pH[4,1] 1.930 0.475 1.178 1.859 3.183
beta1_pH[5,1] 2.357 0.483 1.551 2.315 3.481
beta1_pH[6,1] 2.295 0.898 1.156 2.045 4.493
beta1_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,1] 2.814 1.063 1.284 2.524 4.845
beta1_pH[9,1] 2.178 0.681 1.199 2.044 4.003
beta1_pH[10,1] 2.143 0.329 1.570 2.120 2.856
beta1_pH[11,1] 3.207 0.438 2.393 3.189 4.113
beta1_pH[12,1] 2.880 0.290 2.331 2.885 3.456
beta1_pH[13,1] 3.454 0.245 2.978 3.446 3.955
beta1_pH[14,1] 2.835 0.303 2.252 2.828 3.455
beta1_pH[15,1] 2.407 0.271 1.898 2.400 2.953
beta1_pH[16,1] 4.029 0.290 3.497 4.017 4.632
beta1_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[1,3] 3.457 1.670 1.214 3.008 7.467
beta1_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,3] 4.905 2.361 2.150 4.132 11.485
beta1_pH[9,3] 8.411 2.646 4.119 8.050 14.187
beta1_pH[10,3] 7.330 2.811 3.165 6.822 13.480
beta1_pH[11,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[11,4] 1.920 0.710 0.439 1.909 3.330
beta1_pH[12,4] 3.878 1.091 1.815 3.842 6.374
beta1_pH[13,4] 2.417 0.826 1.004 2.338 4.344
beta1_pH[14,4] 2.677 0.662 1.532 2.629 4.170
beta1_pH[15,4] 1.063 0.754 0.115 0.977 2.645
beta1_pH[16,4] 3.735 0.872 2.195 3.663 5.715
beta1_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,1] 1.277 0.880 0.317 0.957 3.408
beta2_pH[2,1] 1.629 0.965 0.325 1.444 3.423
beta2_pH[3,1] 1.388 0.855 0.210 1.253 2.954
beta2_pH[4,1] 1.143 0.809 0.209 0.866 2.924
beta2_pH[5,1] 2.095 1.751 0.282 1.640 7.001
beta2_pH[6,1] 0.905 1.453 0.105 0.369 5.411
beta2_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,1] 0.573 1.067 0.104 0.278 3.388
beta2_pH[9,1] 0.861 1.255 0.139 0.468 4.864
beta2_pH[10,1] 1.250 1.307 0.310 0.819 5.330
beta2_pH[11,1] 0.477 0.306 0.206 0.391 1.237
beta2_pH[12,1] 1.562 0.813 0.614 1.301 3.664
beta2_pH[13,1] 1.010 0.413 0.550 0.913 2.123
beta2_pH[14,1] 1.154 0.704 0.457 0.935 3.422
beta2_pH[15,1] 1.033 0.649 0.367 0.845 3.080
beta2_pH[16,1] 1.063 0.331 0.645 1.002 1.860
beta2_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,3] 2.182 3.867 0.109 0.502 14.492
beta2_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,3] 2.405 2.018 0.570 1.534 8.050
beta2_pH[9,3] 0.185 0.125 0.080 0.159 0.454
beta2_pH[10,3] 0.302 0.352 0.099 0.225 0.931
beta2_pH[11,3] -8.289 6.039 -22.911 -7.075 -0.388
beta2_pH[12,3] -8.520 6.065 -22.945 -7.306 -0.458
beta2_pH[13,3] -8.395 6.047 -22.721 -7.127 -0.450
beta2_pH[14,3] -8.377 6.011 -22.500 -7.092 -0.375
beta2_pH[15,3] -8.414 6.057 -22.816 -7.182 -0.357
beta2_pH[16,3] -8.389 6.032 -22.820 -7.164 -0.426
beta2_pH[11,4] -5.097 4.179 -16.487 -3.821 -0.696
beta2_pH[12,4] -5.655 4.202 -16.728 -4.485 -0.950
beta2_pH[13,4] -4.943 4.235 -17.160 -3.446 -0.738
beta2_pH[14,4] -4.958 4.041 -16.033 -3.552 -0.864
beta2_pH[15,4] -5.425 4.243 -16.946 -4.132 -0.856
beta2_pH[16,4] -5.411 3.957 -16.358 -4.110 -1.318
beta2_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta3_pH[1,1] 34.982 1.123 32.872 35.003 37.377
beta3_pH[2,1] 32.543 1.324 30.144 32.493 35.355
beta3_pH[3,1] 33.052 1.418 30.467 32.990 36.288
beta3_pH[4,1] 32.621 1.930 29.686 32.376 37.245
beta3_pH[5,1] 27.537 1.081 25.857 27.445 30.339
beta3_pH[6,1] 33.785 3.435 27.030 33.396 41.485
beta3_pH[7,1] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,1] 36.142 3.846 28.749 36.429 42.741
beta3_pH[9,1] 32.652 3.020 27.725 32.132 39.806
beta3_pH[10,1] 32.470 1.160 30.252 32.479 34.787
beta3_pH[11,1] 31.259 1.086 29.256 31.285 33.360
beta3_pH[12,1] 30.814 0.472 29.865 30.828 31.712
beta3_pH[13,1] 31.411 0.516 30.457 31.422 32.463
beta3_pH[14,1] 32.492 0.770 31.157 32.449 34.094
beta3_pH[15,1] 30.896 0.824 29.292 30.915 32.487
beta3_pH[16,1] 30.129 0.342 29.453 30.125 30.800
beta3_pH[1,2] 34.667 5.433 25.449 34.875 43.445
beta3_pH[2,2] 34.218 5.416 25.464 34.168 43.495
beta3_pH[3,2] 30.525 7.237 18.726 30.527 42.452
beta3_pH[4,2] 30.489 7.182 18.753 30.296 42.493
beta3_pH[5,2] 30.606 7.212 18.604 30.769 42.365
beta3_pH[6,2] 30.290 7.160 18.794 30.173 42.166
beta3_pH[7,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[9,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[10,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[11,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[12,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[13,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[14,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[15,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[16,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[1,3] 39.804 2.563 33.328 40.153 43.722
beta3_pH[2,3] 34.499 5.443 25.541 34.406 43.439
beta3_pH[3,3] 30.600 7.126 18.630 30.495 42.389
beta3_pH[4,3] 30.433 7.254 18.531 30.234 42.332
beta3_pH[5,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[6,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[7,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,3] 42.144 0.847 40.988 41.885 43.878
beta3_pH[9,3] 35.897 1.609 34.525 35.282 40.565
beta3_pH[10,3] 37.862 2.568 34.628 37.135 43.445
beta3_pH[11,3] 38.465 2.592 34.171 38.479 42.737
beta3_pH[12,3] 38.546 2.563 34.259 38.606 42.777
beta3_pH[13,3] 38.431 2.629 34.160 38.390 42.783
beta3_pH[14,3] 38.453 2.558 34.218 38.423 42.688
beta3_pH[15,3] 38.521 2.610 34.228 38.524 42.817
beta3_pH[16,3] 38.463 2.604 34.247 38.477 42.790
beta3_pH[11,4] 40.501 1.192 38.200 40.417 43.082
beta3_pH[12,4] 42.185 0.730 40.348 42.371 43.027
beta3_pH[13,4] 40.815 0.986 38.985 40.823 42.915
beta3_pH[14,4] 40.174 0.711 38.836 40.136 41.587
beta3_pH[15,4] 41.475 1.490 37.677 41.636 43.817
beta3_pH[16,4] 41.684 0.461 40.800 41.726 42.581
beta3_pH[11,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[12,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[13,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[14,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[15,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[16,5] 25.000 0.000 25.000 25.000 25.000
beta4_pH[1,1] 0.263 0.021 0.223 0.263 0.304
beta4_pH[2,1] 0.211 0.005 0.201 0.211 0.221
beta4_pH[3,1] 0.172 0.007 0.158 0.172 0.187
beta4_pH[4,1] 0.294 0.011 0.271 0.294 0.316
beta4_pH[5,1] 0.412 0.264 0.083 0.374 0.934
beta4_pH[6,1] 0.764 0.332 0.315 0.726 1.405
beta4_pH[7,1] 0.842 0.329 0.356 0.791 1.615
beta4_pH[8,1] 0.157 0.069 0.041 0.150 0.308
beta4_pH[9,1] 0.215 0.086 0.065 0.210 0.397
beta4_pH[10,1] 0.001 0.001 0.000 0.001 0.003
beta4_pH[11,1] 0.446 0.012 0.424 0.446 0.468
beta4_pH[12,1] 0.643 0.023 0.600 0.643 0.691
beta4_pH[13,1] 0.630 0.014 0.603 0.630 0.657
beta4_pH[14,1] 0.670 0.031 0.609 0.669 0.733
beta4_pH[15,1] 0.797 0.021 0.757 0.796 0.838
beta4_pH[16,1] 1.214 0.041 1.133 1.214 1.292
beta4_pH[1,2] 0.298 0.096 0.128 0.294 0.498
beta4_pH[2,2] 0.514 0.023 0.471 0.513 0.559
beta4_pH[3,2] 0.446 0.012 0.424 0.446 0.470
beta4_pH[4,2] 0.654 0.037 0.584 0.655 0.730
beta4_pH[5,2] 0.727 0.997 0.099 0.524 2.478
beta4_pH[6,2] 0.916 0.793 0.199 0.746 2.709
beta4_pH[7,2] 0.894 0.444 0.352 0.815 1.933
beta4_pH[8,2] 0.845 0.205 0.474 0.832 1.286
beta4_pH[9,2] 1.158 1.205 0.288 0.871 3.651
beta4_pH[10,2] 0.456 0.051 0.359 0.454 0.557
beta4_pH[11,2] 1.271 0.051 1.174 1.271 1.370
beta4_pH[12,2] 0.609 0.110 0.408 0.604 0.842
beta4_pH[13,2] 0.954 0.045 0.866 0.952 1.045
beta4_pH[14,2] 0.933 0.062 0.813 0.932 1.056
beta4_pH[15,2] 1.009 0.025 0.961 1.009 1.060
beta4_pH[16,2] 1.389 0.075 1.243 1.389 1.537
beta4_pH[1,3] 0.243 0.078 0.104 0.240 0.411
beta4_pH[2,3] 0.050 0.017 0.019 0.050 0.086
beta4_pH[3,3] 0.052 0.010 0.033 0.052 0.071
beta4_pH[4,3] 0.166 0.030 0.107 0.166 0.227
beta4_pH[5,3] 0.744 0.935 0.115 0.547 2.665
beta4_pH[6,3] 0.871 0.694 0.193 0.707 2.561
beta4_pH[7,3] 0.884 0.425 0.350 0.817 1.795
beta4_pH[8,3] 0.394 0.184 0.101 0.374 0.822
beta4_pH[9,3] 0.918 0.846 0.167 0.706 2.833
beta4_pH[10,3] 0.056 0.038 0.003 0.050 0.143
beta4_pH[11,3] 0.621 0.264 0.262 0.572 1.279
beta4_pH[12,3] 0.735 0.254 0.383 0.691 1.341
beta4_pH[13,3] 0.805 0.377 0.315 0.729 1.759
beta4_pH[14,3] 0.862 0.216 0.504 0.837 1.348
beta4_pH[15,3] 0.794 0.201 0.469 0.774 1.242
beta4_pH[16,3] 0.833 0.381 0.317 0.764 1.754
beta4_pH[1,4] 0.639 0.379 0.200 0.554 1.609
beta4_pH[2,4] 0.818 0.929 0.112 0.588 2.994
beta4_pH[3,4] 0.894 0.984 0.110 0.619 3.288
beta4_pH[4,4] 0.746 0.590 0.216 0.627 2.011
beta4_pH[5,4] 0.782 1.003 0.108 0.556 2.893
beta4_pH[6,4] 0.919 0.829 0.217 0.736 2.623
beta4_pH[7,4] 0.893 0.452 0.343 0.817 1.909
beta4_pH[8,4] 0.742 0.715 0.144 0.588 2.230
beta4_pH[9,4] 0.771 0.918 0.119 0.555 2.834
beta4_pH[10,4] 1.422 2.822 0.007 0.378 9.833
beta4_pH[11,4] 0.427 0.019 0.390 0.427 0.464
beta4_pH[12,4] 0.535 0.044 0.448 0.534 0.624
beta4_pH[13,4] 0.836 0.033 0.772 0.835 0.898
beta4_pH[14,4] 0.876 0.061 0.761 0.875 1.000
beta4_pH[15,4] 0.617 0.018 0.582 0.618 0.653
beta4_pH[16,4] 0.664 0.043 0.579 0.663 0.748
beta4_pH[1,5] 0.628 0.349 0.206 0.555 1.516
beta4_pH[2,5] 0.804 0.837 0.107 0.575 3.047
beta4_pH[3,5] 0.916 1.082 0.105 0.607 3.589
beta4_pH[4,5] 0.731 0.479 0.213 0.617 1.947
beta4_pH[5,5] 0.756 0.878 0.098 0.559 2.582
beta4_pH[6,5] 0.925 0.775 0.212 0.746 2.784
beta4_pH[7,5] 0.901 0.422 0.345 0.819 1.891
beta4_pH[8,5] 0.740 0.646 0.140 0.596 2.376
beta4_pH[9,5] 0.805 1.080 0.121 0.565 2.747
beta4_pH[10,5] 1.600 3.339 0.008 0.400 12.205
beta4_pH[11,5] 0.309 0.082 0.159 0.305 0.483
beta4_pH[12,5] 0.873 0.299 0.455 0.813 1.620
beta4_pH[13,5] 2.563 0.358 1.889 2.561 3.316
beta4_pH[14,5] 0.851 0.196 0.529 0.833 1.290
beta4_pH[15,5] 0.931 0.079 0.795 0.924 1.110
beta4_pH[16,5] 0.966 0.391 0.389 0.901 1.938
beta0_pelagic[1] 2.268 0.036 2.197 2.268 2.341
beta0_pelagic[2] 1.538 0.007 1.524 1.538 1.552
beta0_pelagic[3] 0.206 0.261 -0.421 0.273 0.578
beta0_pelagic[4] 0.184 0.232 -0.395 0.207 0.573
beta0_pelagic[5] 1.443 0.157 1.144 1.441 1.761
beta0_pelagic[6] 1.356 0.061 1.246 1.354 1.479
beta0_pelagic[7] 1.471 0.039 1.399 1.470 1.550
beta0_pelagic[8] 1.671 0.043 1.590 1.670 1.754
beta0_pelagic[9] 2.776 0.063 2.655 2.775 2.909
beta0_pelagic[10] 2.603 0.031 2.544 2.602 2.663
beta0_pelagic[11] -0.129 0.197 -0.586 -0.098 0.192
beta0_pelagic[12] 1.800 0.035 1.736 1.799 1.870
beta0_pelagic[13] -1.246 0.633 -2.637 -1.179 -0.202
beta0_pelagic[14] -0.260 0.162 -0.607 -0.252 0.036
beta0_pelagic[15] -0.211 0.035 -0.301 -0.196 -0.180
beta0_pelagic[16] -0.597 0.386 -1.360 -0.557 0.000
beta1_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[3] 0.741 0.285 0.296 0.693 1.454
beta1_pelagic[4] 0.868 0.259 0.424 0.846 1.492
beta1_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[11] 1.141 0.275 0.687 1.101 1.768
beta1_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[13] 3.399 0.802 1.738 3.406 4.884
beta1_pelagic[14] 1.707 0.396 1.099 1.641 2.712
beta1_pelagic[15] 0.102 0.168 0.000 0.017 0.625
beta1_pelagic[16] 1.767 0.611 0.844 1.699 3.005
beta2_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[3] 6.521 6.156 0.215 4.791 22.409
beta2_pelagic[4] 6.985 5.907 0.275 5.546 22.028
beta2_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[11] 1.803 1.546 0.282 1.295 5.938
beta2_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[13] 0.147 0.059 0.078 0.133 0.301
beta2_pelagic[14] 0.924 0.841 0.230 0.650 3.457
beta2_pelagic[15] 2.681 1.921 0.167 2.306 7.188
beta2_pelagic[16] 0.629 0.897 0.112 0.295 3.748
beta3_pelagic[1] 28.826 3.905 22.414 28.888 35.162
beta3_pelagic[2] 28.712 3.937 22.333 28.722 35.158
beta3_pelagic[3] 27.979 2.763 23.213 28.466 32.580
beta3_pelagic[4] 25.940 1.283 23.047 26.214 28.566
beta3_pelagic[5] 29.593 6.629 18.619 29.803 40.509
beta3_pelagic[6] 29.337 6.650 18.546 29.259 40.358
beta3_pelagic[7] 29.536 6.595 18.696 29.467 40.453
beta3_pelagic[8] 29.649 6.677 18.662 29.783 40.437
beta3_pelagic[9] 29.479 6.630 18.659 29.418 40.454
beta3_pelagic[10] 29.472 6.663 18.537 29.453 40.373
beta3_pelagic[11] 30.218 1.319 27.739 30.147 33.229
beta3_pelagic[12] 30.000 0.000 30.000 30.000 30.000
beta3_pelagic[13] 31.779 4.993 23.582 31.552 41.138
beta3_pelagic[14] 34.075 1.596 31.384 33.926 38.032
beta3_pelagic[15] 34.289 5.737 23.508 35.089 41.770
beta3_pelagic[16] 30.717 2.865 25.648 30.541 37.691
beta0_yellow[1] -0.485 0.080 -0.637 -0.485 -0.335
beta0_yellow[2] 0.537 0.090 0.316 0.547 0.683
beta0_yellow[3] -0.208 0.125 -0.457 -0.214 0.014
beta0_yellow[4] 0.874 0.274 0.105 0.952 1.194
beta0_yellow[5] 0.248 0.333 -0.471 0.276 0.818
beta0_yellow[6] 0.817 0.104 0.620 0.817 1.020
beta0_yellow[7] 1.013 0.101 0.820 1.010 1.214
beta0_yellow[8] 1.017 0.093 0.829 1.017 1.201
beta0_yellow[9] 0.442 0.100 0.248 0.444 0.637
beta0_yellow[10] 0.539 0.071 0.402 0.537 0.677
beta0_yellow[11] -0.191 0.118 -0.454 -0.180 0.017
beta0_yellow[12] -1.093 0.096 -1.279 -1.096 -0.909
beta0_yellow[13] -1.354 0.173 -1.723 -1.343 -1.057
beta0_yellow[14] -0.259 0.167 -0.611 -0.237 -0.015
beta0_yellow[15] -1.235 0.121 -1.509 -1.224 -1.029
beta0_yellow[16] -0.255 0.094 -0.437 -0.256 -0.066
beta1_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[2] 1.012 0.242 0.577 1.004 1.496
beta1_yellow[3] 0.469 0.248 0.019 0.482 0.917
beta1_yellow[4] 1.384 0.639 0.657 1.206 3.358
beta1_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[11] 0.225 0.299 0.001 0.057 0.987
beta1_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[13] 0.917 0.892 0.023 0.717 3.865
beta1_yellow[14] 0.643 0.931 0.012 0.369 3.709
beta1_yellow[15] 0.153 0.199 0.011 0.073 0.670
beta1_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[2] -7.649 5.768 -21.592 -6.366 -0.375
beta2_yellow[3] -7.833 5.704 -21.282 -6.722 -0.438
beta2_yellow[4] -5.458 5.878 -19.934 -3.335 -0.116
beta2_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[11] -7.772 5.737 -21.002 -6.615 -0.379
beta2_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[13] -6.247 5.370 -19.430 -4.723 -0.315
beta2_yellow[14] -6.888 5.194 -18.902 -6.102 -0.266
beta2_yellow[15] -7.413 5.514 -20.648 -6.285 -0.366
beta2_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta3_yellow[1] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[2] 28.982 1.245 27.324 28.704 32.531
beta3_yellow[3] 32.404 3.317 23.242 32.653 38.661
beta3_yellow[4] 29.801 3.401 23.073 28.637 35.903
beta3_yellow[5] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[6] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[7] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[8] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[9] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[10] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[11] 35.251 3.785 30.352 34.210 43.378
beta3_yellow[12] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[13] 34.066 2.712 29.476 34.052 40.664
beta3_yellow[14] 34.236 3.690 29.190 33.761 43.085
beta3_yellow[15] 36.436 3.837 30.279 36.295 43.662
beta3_yellow[16] 35.000 0.000 35.000 35.000 35.000
beta0_black[1] 0.013 0.074 -0.130 0.012 0.161
beta0_black[2] 1.855 0.032 1.792 1.855 1.919
beta0_black[3] 1.240 0.038 1.164 1.240 1.314
beta0_black[4] 2.425 0.047 2.333 2.426 2.512
beta0_black[5] 5.194 2.466 1.829 4.677 11.488
beta0_black[6] 5.270 2.229 2.663 4.659 11.396
beta0_black[7] 4.100 2.123 1.591 3.537 10.020
beta0_black[8] 0.885 0.216 0.481 0.879 1.345
beta0_black[9] 2.807 0.265 2.306 2.807 3.335
beta0_black[10] 1.288 0.050 1.193 1.286 1.386
beta0_black[11] 2.978 0.052 2.876 2.978 3.079
beta0_black[12] 4.916 0.194 4.538 4.919 5.297
beta0_black[13] 0.445 0.332 0.013 0.369 1.226
beta0_black[14] 2.852 0.105 2.645 2.853 3.059
beta0_black[15] 0.722 0.052 0.620 0.723 0.824
beta0_black[16] 3.773 0.074 3.630 3.773 3.919
beta2_black[1] 0.000 0.000 0.000 0.000 0.000
beta2_black[2] 0.000 0.000 0.000 0.000 0.000
beta2_black[3] 0.000 0.000 0.000 0.000 0.000
beta2_black[4] 0.000 0.000 0.000 0.000 0.000
beta2_black[5] 0.000 0.000 0.000 0.000 0.000
beta2_black[6] 0.000 0.000 0.000 0.000 0.000
beta2_black[7] 0.000 0.000 0.000 0.000 0.000
beta2_black[8] 0.000 0.000 0.000 0.000 0.000
beta2_black[9] 0.000 0.000 0.000 0.000 0.000
beta2_black[10] 0.000 0.000 0.000 0.000 0.000
beta2_black[11] 0.000 0.000 0.000 0.000 0.000
beta2_black[12] 0.000 0.000 0.000 0.000 0.000
beta2_black[13] -2.573 1.806 -7.002 -2.173 -0.281
beta2_black[14] 0.000 0.000 0.000 0.000 0.000
beta2_black[15] 0.000 0.000 0.000 0.000 0.000
beta2_black[16] 0.000 0.000 0.000 0.000 0.000
beta3_black[1] 25.000 0.000 25.000 25.000 25.000
beta3_black[2] 25.000 0.000 25.000 25.000 25.000
beta3_black[3] 25.000 0.000 25.000 25.000 25.000
beta3_black[4] 25.000 0.000 25.000 25.000 25.000
beta3_black[5] 25.000 0.000 25.000 25.000 25.000
beta3_black[6] 25.000 0.000 25.000 25.000 25.000
beta3_black[7] 25.000 0.000 25.000 25.000 25.000
beta3_black[8] 25.000 0.000 25.000 25.000 25.000
beta3_black[9] 25.000 0.000 25.000 25.000 25.000
beta3_black[10] 25.000 0.000 25.000 25.000 25.000
beta3_black[11] 25.000 0.000 25.000 25.000 25.000
beta3_black[12] 25.000 0.000 25.000 25.000 25.000
beta3_black[13] 37.450 3.117 28.536 38.115 41.077
beta3_black[14] 25.000 0.000 25.000 25.000 25.000
beta3_black[15] 25.000 0.000 25.000 25.000 25.000
beta3_black[16] 25.000 0.000 25.000 25.000 25.000
beta4_black[1] -0.327 0.091 -0.509 -0.326 -0.146
beta4_black[2] 0.258 0.053 0.161 0.258 0.361
beta4_black[3] -0.721 0.064 -0.843 -0.721 -0.597
beta4_black[4] 0.627 0.130 0.376 0.623 0.879
beta4_black[5] 0.800 1.616 -1.185 0.506 4.565
beta4_black[6] 0.866 1.582 -0.868 0.557 4.159
beta4_black[7] 0.791 1.550 -0.977 0.521 4.136
beta4_black[8] 0.084 0.282 -0.530 0.095 0.605
beta4_black[9] 1.039 0.798 -0.063 0.888 3.032
beta4_black[10] 0.075 0.066 -0.055 0.074 0.204
beta4_black[11] -0.608 0.065 -0.731 -0.610 -0.480
beta4_black[12] 1.235 0.288 0.690 1.232 1.799
beta4_black[13] -1.003 0.078 -1.156 -1.004 -0.852
beta4_black[14] 0.217 0.181 -0.136 0.219 0.566
beta4_black[15] 0.002 0.046 -0.086 0.002 0.093
beta4_black[16] -0.414 0.116 -0.641 -0.413 -0.187
beta0_dsr[11] -0.338 0.009 -0.357 -0.338 -0.320
beta0_dsr[12] 1.420 0.044 1.335 1.421 1.507
beta0_dsr[13] 0.143 0.021 0.102 0.143 0.185
beta0_dsr[14] -0.128 0.017 -0.160 -0.128 -0.095
beta0_dsr[15] 0.196 0.014 0.168 0.196 0.224
beta0_dsr[16] 0.208 0.015 0.179 0.208 0.239
beta1_dsr[11] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[12] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[13] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[14] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[15] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[16] 0.000 0.000 0.000 0.000 0.000
beta4_dsr[11] 0.548 0.036 0.479 0.549 0.618
beta4_dsr[12] 0.788 0.232 0.348 0.786 1.256
beta4_dsr[13] -0.078 0.040 -0.159 -0.079 0.002
beta4_dsr[14] -0.214 0.073 -0.360 -0.213 -0.075
beta4_dsr[15] 0.227 0.022 0.184 0.226 0.268
beta4_dsr[16] -0.098 0.045 -0.186 -0.098 -0.011
beta0_slope[11] -2.198 0.021 -2.240 -2.198 -2.158
beta0_slope[12] -5.150 0.176 -5.480 -5.152 -4.780
beta0_slope[13] -1.473 0.035 -1.542 -1.473 -1.407
beta0_slope[14] -2.685 0.048 -2.780 -2.684 -2.595
beta0_slope[15] -1.590 0.023 -1.636 -1.590 -1.543
beta0_slope[16] -3.015 0.047 -3.109 -3.015 -2.925
beta1_slope[11] 0.000 0.000 0.000 0.000 0.000
beta1_slope[12] 0.000 0.000 0.000 0.000 0.000
beta1_slope[13] 0.000 0.000 0.000 0.000 0.000
beta1_slope[14] 0.000 0.000 0.000 0.000 0.000
beta1_slope[15] 0.000 0.000 0.000 0.000 0.000
beta1_slope[16] 0.000 0.000 0.000 0.000 0.000
beta4_slope[11] -0.345 0.075 -0.495 -0.344 -0.204
beta4_slope[12] -4.274 0.825 -6.043 -4.211 -2.843
beta4_slope[13] 0.309 0.053 0.206 0.307 0.418
beta4_slope[14] -0.188 0.147 -0.474 -0.186 0.098
beta4_slope[15] -0.364 0.034 -0.431 -0.364 -0.299
beta4_slope[16] 0.222 0.105 0.014 0.224 0.429
beta5_slope[11] 3.071 0.108 2.868 3.070 3.287
beta5_slope[12] 1.576 0.417 0.845 1.558 2.436
beta5_slope[13] 0.847 0.079 0.690 0.847 1.000
beta5_slope[14] 3.781 0.369 3.130 3.768 4.550
beta5_slope[15] 1.415 0.072 1.276 1.414 1.561
beta5_slope[16] 2.681 0.227 2.265 2.672 3.157
tau_comp[1] 3.790 0.287 3.252 3.782 4.377
tau_comp[2] 2.859 0.261 2.394 2.848 3.404
tau_comp[3] 0.872 0.076 0.731 0.868 1.029
tau_comp[4] 1.354 0.213 0.990 1.335 1.824
tau_comp[5] 1.206 0.172 0.906 1.196 1.570
sigma_H[1] 0.204 0.047 0.116 0.201 0.303
sigma_H[2] 0.174 0.028 0.127 0.171 0.237
sigma_H[3] 0.189 0.040 0.114 0.186 0.277
sigma_H[4] 0.463 0.079 0.331 0.454 0.638
sigma_H[5] 0.966 0.202 0.598 0.955 1.387
sigma_H[6] 0.541 0.195 0.160 0.535 0.933
sigma_H[7] 0.310 0.061 0.215 0.301 0.450
sigma_H[8] 0.417 0.089 0.262 0.410 0.603
sigma_H[9] 0.513 0.117 0.329 0.499 0.779
sigma_H[10] 0.209 0.039 0.144 0.205 0.294
sigma_H[11] 0.266 0.045 0.190 0.261 0.368
sigma_H[12] 0.382 0.148 0.202 0.333 0.718
sigma_H[13] 0.221 0.036 0.159 0.218 0.300
sigma_H[14] 0.499 0.089 0.346 0.492 0.699
sigma_H[15] 0.231 0.038 0.167 0.227 0.319
sigma_H[16] 0.213 0.040 0.144 0.210 0.306
lambda_H[1] 2.146 3.501 0.078 0.962 11.581
lambda_H[2] 5.298 5.184 0.495 3.840 18.472
lambda_H[3] 5.119 7.038 0.288 2.698 24.001
lambda_H[4] 0.006 0.004 0.001 0.005 0.018
lambda_H[5] 2.637 6.432 0.031 0.670 19.704
lambda_H[6] 6.407 13.241 0.008 0.412 45.591
lambda_H[7] 0.012 0.008 0.002 0.010 0.033
lambda_H[8] 7.857 9.894 0.092 4.261 35.938
lambda_H[9] 0.015 0.010 0.003 0.013 0.040
lambda_H[10] 0.223 0.278 0.025 0.152 0.874
lambda_H[11] 0.225 0.332 0.010 0.110 1.070
lambda_H[12] 3.610 5.223 0.175 2.068 15.786
lambda_H[13] 3.241 3.154 0.257 2.343 11.441
lambda_H[14] 3.456 4.108 0.239 2.182 15.337
lambda_H[15] 0.018 0.019 0.003 0.014 0.059
lambda_H[16] 0.534 0.771 0.037 0.282 2.405
mu_lambda_H[1] 4.086 1.898 1.026 3.886 8.294
mu_lambda_H[2] 3.657 1.966 0.465 3.491 7.793
mu_lambda_H[3] 3.414 1.869 0.705 3.120 7.707
sigma_lambda_H[1] 8.057 4.252 1.698 7.493 17.679
sigma_lambda_H[2] 8.036 4.702 0.741 7.440 18.199
sigma_lambda_H[3] 6.288 4.095 0.956 5.435 16.440
beta_H[1,1] 6.475 1.341 3.310 6.688 8.421
beta_H[2,1] 9.778 0.596 8.506 9.813 10.854
beta_H[3,1] 7.815 0.787 5.879 7.898 9.136
beta_H[4,1] 10.471 7.467 -4.492 10.588 25.406
beta_H[5,1] -0.285 2.520 -5.361 -0.158 4.055
beta_H[6,1] 3.175 4.330 -7.615 4.769 8.239
beta_H[7,1] 1.106 5.618 -10.500 1.370 11.519
beta_H[8,1] 1.287 3.878 -2.681 1.228 3.506
beta_H[9,1] 13.668 5.789 2.409 13.565 25.340
beta_H[10,1] 6.831 1.850 3.095 6.932 10.514
beta_H[11,1] 4.868 3.512 -2.728 5.469 9.851
beta_H[12,1] 2.626 1.089 0.677 2.549 5.054
beta_H[13,1] 9.073 0.851 7.243 9.127 10.522
beta_H[14,1] 2.153 1.016 0.213 2.141 4.197
beta_H[15,1] -7.050 3.310 -13.123 -7.232 0.003
beta_H[16,1] 3.793 2.845 -0.831 3.440 10.332
beta_H[1,2] 7.780 0.270 7.217 7.783 8.292
beta_H[2,2] 9.885 0.161 9.570 9.886 10.206
beta_H[3,2] 8.792 0.200 8.391 8.795 9.191
beta_H[4,2] 3.058 1.433 0.342 3.042 5.926
beta_H[5,2] 1.757 0.980 -0.207 1.754 3.646
beta_H[6,2] 6.076 1.074 3.638 6.208 7.861
beta_H[7,2] 2.418 1.058 0.515 2.382 4.633
beta_H[8,2] 2.998 1.116 1.222 3.154 4.243
beta_H[9,2] 3.184 1.122 1.048 3.187 5.535
beta_H[10,2] 8.038 0.368 7.283 8.052 8.740
beta_H[11,2] 9.823 0.633 8.861 9.715 11.200
beta_H[12,2] 3.935 0.370 3.245 3.928 4.701
beta_H[13,2] 9.143 0.245 8.687 9.140 9.637
beta_H[14,2] 4.012 0.340 3.349 4.008 4.689
beta_H[15,2] 11.539 0.597 10.290 11.562 12.659
beta_H[16,2] 4.476 0.819 2.974 4.449 6.140
beta_H[1,3] 8.448 0.263 7.960 8.440 8.968
beta_H[2,3] 9.859 0.147 9.566 9.861 10.157
beta_H[3,3] 9.459 0.165 9.147 9.452 9.809
beta_H[4,3] -2.165 0.876 -3.919 -2.175 -0.452
beta_H[5,3] 3.891 0.602 2.698 3.894 5.097
beta_H[6,3] 8.347 1.181 6.679 8.002 10.904
beta_H[7,3] -2.573 0.621 -3.810 -2.569 -1.387
beta_H[8,3] 5.303 0.527 4.701 5.227 6.440
beta_H[9,3] -2.723 0.704 -4.174 -2.703 -1.408
beta_H[10,3] 8.542 0.280 7.975 8.541 9.095
beta_H[11,3] 8.505 0.287 7.887 8.528 9.005
beta_H[12,3] 5.258 0.319 4.518 5.296 5.787
beta_H[13,3] 8.861 0.175 8.499 8.867 9.195
beta_H[14,3] 5.771 0.270 5.173 5.794 6.251
beta_H[15,3] 10.266 0.294 9.705 10.260 10.867
beta_H[16,3] 6.000 0.576 4.831 6.028 7.036
beta_H[1,4] 8.089 0.185 7.712 8.099 8.426
beta_H[2,4] 9.854 0.131 9.581 9.861 10.092
beta_H[3,4] 9.949 0.150 9.624 9.956 10.221
beta_H[4,4] 11.802 0.473 10.866 11.821 12.711
beta_H[5,4] 5.577 0.728 4.368 5.503 7.203
beta_H[6,4] 7.234 0.912 5.314 7.456 8.514
beta_H[7,4] 8.545 0.343 7.887 8.542 9.217
beta_H[8,4] 6.724 0.267 6.216 6.736 7.170
beta_H[9,4] 7.434 0.443 6.580 7.429 8.329
beta_H[10,4] 7.529 0.228 7.101 7.524 8.007
beta_H[11,4] 9.505 0.190 9.130 9.503 9.885
beta_H[12,4] 7.251 0.198 6.860 7.246 7.670
beta_H[13,4] 9.055 0.143 8.763 9.059 9.331
beta_H[14,4] 7.818 0.211 7.419 7.809 8.260
beta_H[15,4] 9.467 0.220 9.026 9.467 9.917
beta_H[16,4] 9.417 0.235 8.969 9.410 9.887
beta_H[1,5] 8.898 0.151 8.590 8.899 9.192
beta_H[2,5] 10.676 0.100 10.496 10.672 10.887
beta_H[3,5] 10.764 0.159 10.467 10.759 11.081
beta_H[4,5] 8.348 0.491 7.434 8.334 9.369
beta_H[5,5] 5.203 0.592 3.742 5.266 6.185
beta_H[6,5] 8.958 0.629 8.012 8.820 10.360
beta_H[7,5] 6.567 0.340 5.892 6.563 7.260
beta_H[8,5] 8.247 0.216 7.902 8.228 8.664
beta_H[9,5] 8.191 0.446 7.269 8.196 9.041
beta_H[10,5] 10.044 0.231 9.583 10.046 10.486
beta_H[11,5] 11.592 0.222 11.156 11.592 12.017
beta_H[12,5] 8.618 0.187 8.248 8.618 8.989
beta_H[13,5] 10.125 0.132 9.882 10.119 10.404
beta_H[14,5] 9.240 0.218 8.840 9.230 9.693
beta_H[15,5] 11.290 0.235 10.834 11.288 11.741
beta_H[16,5] 9.868 0.184 9.480 9.873 10.218
beta_H[1,6] 10.126 0.237 9.744 10.095 10.710
beta_H[2,6] 11.466 0.112 11.246 11.467 11.679
beta_H[3,6] 10.616 0.155 10.279 10.622 10.891
beta_H[4,6] 12.919 0.868 11.171 12.941 14.543
beta_H[5,6] 5.664 0.632 4.549 5.641 6.932
beta_H[6,6] 8.757 0.569 7.552 8.788 9.795
beta_H[7,6] 10.545 0.581 9.404 10.539 11.661
beta_H[8,6] 9.521 0.283 9.057 9.542 9.944
beta_H[9,6] 8.302 0.751 6.907 8.295 9.888
beta_H[10,6] 9.140 0.347 8.422 9.163 9.761
beta_H[11,6] 10.693 0.358 9.941 10.724 11.335
beta_H[12,6] 9.374 0.240 8.908 9.362 9.896
beta_H[13,6] 11.068 0.150 10.769 11.067 11.369
beta_H[14,6] 9.905 0.281 9.348 9.908 10.451
beta_H[15,6] 10.530 0.416 9.710 10.525 11.335
beta_H[16,6] 10.585 0.253 10.062 10.591 11.088
beta_H[1,7] 10.084 1.280 6.771 10.397 11.821
beta_H[2,7] 12.139 0.498 11.057 12.168 13.036
beta_H[3,7] 10.290 0.671 8.786 10.355 11.419
beta_H[4,7] 1.429 4.446 -7.247 1.321 10.297
beta_H[5,7] 5.873 2.051 1.964 5.883 9.801
beta_H[6,7] 7.206 3.169 -0.850 8.394 10.796
beta_H[7,7] 7.665 2.955 1.922 7.691 13.606
beta_H[8,7] 10.874 1.059 9.101 10.851 12.383
beta_H[9,7] 4.893 3.900 -3.128 4.939 12.304
beta_H[10,7] 10.077 1.608 7.318 9.943 13.585
beta_H[11,7] 11.370 1.823 8.120 11.232 15.407
beta_H[12,7] 9.776 0.951 7.609 9.848 11.426
beta_H[13,7] 11.671 0.695 10.072 11.736 12.871
beta_H[14,7] 10.627 0.911 8.724 10.657 12.348
beta_H[15,7] 14.061 2.180 9.774 14.041 18.459
beta_H[16,7] 12.854 1.412 10.482 12.712 16.007
beta0_H[1] 7.945 19.205 -32.424 8.575 45.110
beta0_H[2] 10.403 8.083 -5.977 10.479 26.854
beta0_H[3] 9.666 9.733 -9.892 9.572 30.507
beta0_H[4] 6.721 189.191 -381.789 7.837 409.084
beta0_H[5] 3.657 27.886 -50.864 4.342 58.493
beta0_H[6] 7.297 49.562 -109.532 7.972 119.358
beta0_H[7] 7.757 135.577 -259.646 6.054 287.053
beta0_H[8] 6.501 38.526 -20.023 6.277 29.978
beta0_H[9] 5.125 121.564 -239.506 6.711 244.339
beta0_H[10] 8.803 37.912 -64.924 8.377 85.962
beta0_H[11] 7.811 53.368 -105.054 8.467 113.925
beta0_H[12] 6.519 12.836 -18.381 6.674 32.276
beta0_H[13] 9.542 10.215 -11.168 9.616 29.426
beta0_H[14] 7.152 11.335 -15.960 7.118 29.486
beta0_H[15] 7.152 117.810 -233.170 8.143 235.241
beta0_H[16] 9.400 30.607 -52.358 8.480 75.661